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1.  Introduction 


The  primary  effect  of  cyanide  (CN)  on  cardiac  tissue  is  its  interference  with  the  cellular  energy 
balance,  that  is,  its  inhibition  of  cytochrome  oxidase,  a  mitochondrial  enzyme,  preventing  the 
transfer  of  electrons  from  cytochrome  c  to  molecular  oxygen.  The  impacted  aerobic  metabolism 
results  in  lactic  acid  acidosis  and  reduced  adenosine  triphosphate  (ATP)  concentration,  which 
lead  to  anoxia.  With  a  reduction  in  cellular  energy,  the  following  effects  develop  in  tandem: 
changes  in  calcium  ion  concentration  in  several  cellular  compartments,  calcium  current 
modulation,  and  inotropic  effects  (Way,  1984).  Simultaneously,  a  gradual  decrease  in  action 
potential  duration  is  also  noted.  The  subsidiary  cellular  processes  are  still  being  elucidated,  but 
the  sine  qua  non  for  developing  effective  therapeutic  regimes  is  management  of  these  effects. 

An  in  silico  approach  that  shows  how  available  data  and  models  can  be  used  to  develop  new 
therapies  is  presented  in  this  report. 

Toward  these  ends,  an  electrocardiogram  (ECG)  of  an  individual  under  CN  intoxication  is  taken 
as  the  starting  point  (Wexler  et  ah,  1947).  The  deviation  from  homeostasis,  that  is,  from  normal 
conditions,  is  interpreted  in  terms  of  the  modulation  of  the  membrane  ionic  currents  that 
compose  the  action  potential.  The  nature  of  the  deviation  suggests  the  phannacological 
intervention  required  to  stop  and  possibly  reverse  the  toxicology. 

The  need  for  a  therapy  to  counteract  CN  poisoning  is  especially  urgent  in  situations  where  a 
dosage  of  one  LD50  of  CN  is  encountered,  that  is,  where  the  dosage  of  the  toxin  is  >0.2  mg/kg  of 
body  weight.  In  these  cases,  as  shown  by  Ansell  and  Lewis  (1970),  Baskin  et  al.  (1987),  and 
Baskin  (1991),  the  poison  primarily  affects  the  heart,  expressed  through  its  effect  on  the 
ventricles.  Reminiscent  of  severe  ischemia,  lesions  in  the  heart  tissue  were  noted  by  Suzuki 
(1968).  Upon  exposure  to  the  poison,  the  heart  initially  slows  down  into  the  bradycardia  regime, 
then  in  time  transits  to  arrhythmia.  The  ventricles  are  severely  affected,  expressed  by  symptoms 
of  repolarization  disorders  culminating  in  Torsade  de  Pointes  (TdP)  and  ventricular  fibrillation 
(VF). 

An  overview  of  cardiac  cell  energy  metabolism  is  addressed  in  section  2,  setting  the  stage  for  a 
quantification  of  the  changes  brought  about  by  CN  perturbation  of  cardiac  tissue  metabolism. 
Section  3  delves  into  the  issue  of  metabolic  networks  and  sketches  the  opportunities,  in 
conjunction  with  high  perfonnance  computing  requirements,  to  deepen  our  understanding  of  the 
cellular  processes.  This  is  followed  by  an  interpretation  in  section  4  of  the  salient  changes  in  the 
ECG  of  a  subject  under  severe  CN  intoxication.  Section  5  shows  how  in  silico  techniques  can 
shed  light  on  the  action  potential  modulations  caused  by  CN  and  point  to  the  required 
characteristics  of  an  antidote.  Finally,  this  report  concludes  with  a  discussion  of  the  results  in 
section  6. 
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2.  Cardiac  Cell  Energy  Metabolism 


The  cardiac  myocyte  derives  its  energy  primarily  from  fatty  acids,  which  make  up  to  90%  of  the 
total  energy  required.  The  energy  is  derived  from  glucose  and  lactate  metabolism  (Opie,  1998). 
The  mitochondrion  is  the  source  of  the  ATP  produced  from  adenosine  diphosphate  (ADP)  and 
inorganic  phosphate.  Glucose  oxidation,  though  only  a  minor  part  of  the  total  energy  budget  of 
the  cell,  plays  a  crucial  role  in  ion  homeostasis.  An  elaborate  system  of  enzymes  regulates  the 
glycolytic  pathway.  Glycolysis  converts  glucose  to  pyruvate,  which  in  turn  transforms  to 
acetyl-CoA  and  enters  the  Krebs  cycle.  In  a  diseased  state,  the  breakdown  of  glucose  or 
glycogen  to  lactate  replaces  the  oxidative  metabolism  for  anaerobic  energy  production.  In 
anoxia,  ATP  levels  decrease  with  the  cessation  of  oxidative  metabolism,  and  glycolysis  provides 
limited  amounts  of  energy. 

The  three  stages  of  energy  metabolism  are  (1)  Krebs  cycle  input,  (2)  Krebs  cycle  production  of 
nicotinamide  adenine  dinucleotide  (NADH),  and  (3)  oxidative  phosphorylation  (OxPhos). 
Oxidization  takes  place  by  the  electron  transport  chain,  passing  electrons  to  molecular  oxygen  by 
means  of  cytochrome  oxidase  to  fonn  water. 

The  Krebs  cycle  and  glycolysis  act  as  depositories  for  the  energy  obtained  from  substrate 
catabolism.  The  Krebs  cycle  products  are  oxidized  by  the  electron  transport  chain.  ATP 
production  is  controlled  at  the  substrate  level  by  FiFo-ATPase  and  by  ADP  through  a  negative 
feedback  mechanism.  That  is,  an  increase  in  ATP/ ADP  ratio  inhibits  OxPhos.  Also, 
homeostasis  requires  control  of  magnesium  and  calcium  ion  levels  (Wang  et  ah,  2003). 

Energy  metabolism  is  subject  to  a  number  of  pathologies.  Oxygen  uptake  reduction  caused  by 
CN  intoxication  leads  to  a  slowing  of  OxPhos  and  an  increase  in  NADH  levels.  Also,  cell 
acidification,  increased  ADP  and  creatine  production,  and  decreased  creatine  phosphate 
production  are  observed.  Mitochondrial  myopathies  and  cytosolic  pathologies  have  been 
studied,  shedding  light  on  respiratory  chain  defects  in  electron  transport  and  substrate  defects  in 
fatty  acid  utilization.  Cytochrome  oxidase  deficiency,  respiratory  chain  defects,  and  decreases  in 
enzyme  activity  affect  the  cellular  energy  metabolism.* 

Calcium  balance  controls  mitochondrial  respiration.  Mitochondrial  calcium  acts  like  a  buffer  for 
the  cytosol.  Its  importance  is  underlined  by  the  fact  that  excess  calcium  under  cytosolic  calcium 
overload  conditions  limits  the  proton  gradient,  thereby  limiting  the  synthesis  of  ATP. 

This  short  summary  of  energy  metabolism  has  not  touched  upon  many  of  the  subsidiary  issues, 
including  the  kinetics  of  key  regulatory  enzymes,  substrates,  and  products  (Balaban,  2002;  Jafri 
et  ah,  2001).  Rather,  our  aim  is  to  apprise  the  reader  of  the  need  to  consider  the  transformations 


Letellier  et  al.  (1994,  1998)  have  modeled  the  block  of  cytochrome  oxidase  by  cyanide  using  metabolic  control  analysis. 
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and  components  of  metabolic  networks  when  trying  to  understand  the  processes  accompanying 
CN  intoxication.  Respiration  interference  is  not  the  only  process  that  may  play  a  determining 
role. 


3.  Metabolic  Networks:  Work  in  Progress 


Metabolic  and  protein  networks  of  toxin-intoxicated  tissues  contain  a  wealth  of  information 
needed  for  the  development  of  antidotes  and  therapies.  The  networks  are  best  pictured  as  nodes 
and  links.  The  node  may  be  a  substrate  or  a  protein,  and  the  link  is  the  interaction  between  two 
chemical  nodes — a  chemical  transfonnation.  The  interconnected  nodes  form  a  structure. 

The  last  few  years  have  witnessed  an  explosion  in  the  understanding  of  live  and  inanimate 
network  structures  and  the  attempts  to  relate  them  to  performance.  The  research  of  Barabasi  and 
his  group  (Albert  and  Barabasi,  2001;  Barabasi,  2002;  Jeong  et  al.,  2000)  showed  that  a 
network’s  topology  is  the  key  to  its  dynamic  behavior  and  robustness  when  under  attack.  A 
notable  insight  gained  from  tests  on  43  organisms  is  that  their  structures  are  scale-free,  that  is,  no 
single  node  represents  all  of  the  nodes.  In  fact,  the  node  distribution  follows  a  power  law — a  few 
nodes  are  linked  with  many  neighbors,  forming  hubs,  but  the  vast  majority  have  only  a  few 
neighbors.  The  hubs  also  lend  heterogeneity  to  the  topology.  To  test  the  robustness  of  the 
system  to  perturbations,  links  are  removed  at  random.  Most  deletions  do  not  affect  the  overall 
system  function,  but  hub  removal  does.  In  a  biological  network  as  well  as  the  internet,  hub 
removal  is  catastrophic  in  the  sense  that  the  network  is  no  longer  able  to  function  as  intended. 

For  example,  catalytic  enzymes  may  be  removed  from  an  organism  without  adverse  effect  on  its 
overall  performance,  but  if  any  of  the  hubs  are  removed,  the  consequences  could  be  disastrous. 

Questions  about  the  effect  threat  agents  have  on  an  organism  can  be  refonnulated  in  terms  of 
networks.  Identifying  the  affected  link  gives  clues  about  the  cellular  problem  and  possibly 
points  to  alternative  or  compensatory  links  that  could  be  established  to  overcome  the  network 
fault,  a  fault  that  may  adversely  affect  the  organism’s  ability  to  function. 

Further  means  to  understanding  the  mechanisms  involved  in  the  action  of  toxins  include  the 
study  of  metabolic  network  topology,  the  identification  of  hubs,  and  the  study  of  the  kinetics  that 
regulate  the  links  coupled  with  flux  analysis  (Bhalla  and  Iyengar,  1999).  They  shed  light  on  the 
need  for  the  survival  of  certain  node-link  combinations  and,  if  injured,  the  alternate  routes  or 
actions  that  would  enhance  the  likelihood  of  their  survival.  Along  these  lines,  ancillary  work  has 
been  performed  to  determine  the  effect  of  deleting  (or  turning  off)  certain  genes  on  red  blood  cell 
function  and  Escherichia  coli.  For  the  cardiac  myocyte  network,  ATP  is  a  hub.  Its  deletion  or 
isolation  leads  to  apoptosis,  cell  death.  This  is  explored  further  in  section  5. 
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4.  CN-Caused  Changes  in  the  ECG  of  an  Adult  Male 


A  Lead  I  ECG  from  a  NaCN-gas-intoxicated  adult  male,  at  a  lethal  dose  level  of  0.2  mg/kg  of 
body  weight,  is  shown  in  figure  1 .  The  ECG  traces,  interpreted  in  terms  of  myocyte  membrane 
currents,  give  a  time -resolved  record  of  CN's  effect  on  the  heart.  Some  of  the  significant 
anomalies  and  deviations  from  the  norm  in  this  ECG  include  the  following: 

1.  The  initial  heart  rate  (HR)  is  150  beats/min  (BM),  which  changes  to  a  rate  of  37.5  BM  at 
time  t  =  3  min,  signaling  bradycardia.  At  t  =  5  min,  the  rate  increases  to  54.5  BM,  and  at 
t  =  10  min,  before  severe  arrhythmia  develops,  it  increases  to  85.7  BM. 

2.  Morphological  changes  in  the  T-wave  (the  repolarization  of  the  ventricles)  are  apparent 
from  t  =  3  min  on.  There  is  a  notable  widening  of  the  base  of  the  wave  and  a  fusing  of  the 
T-wave  with  the  QRS-complex. 

3.  The  ST-segment  represents  the  early  part  of  the  ventricular  repolarization.  Under  nonnal 
conditions,  its  duration  is  <0.2  s  depending  on  the  HR.  An  ST-segment  changed  by  more 
than  1.0  mm  from  the  horizontal,  homeostatic  position  signals  widespread  myocardial 
injury,  usually  of  the  transmural  kind  (Li  et  ah,  2000;  Yan  and  Antzelevitch,  1999).  Such  a 
change  is  evident  in  figure  1 . 

4.  The  P-wave  is  missing.  By  t  =  3  min,  the  P-wave  has  disappeared,  indicating  no  atrial 
activity.  Absence  of  the  P-wave  persists  from  that  time  on. 

Table  1  summarizes  some  of  the  deviations  from  the  norm  found  in  the  CN-intoxicated 
individual's  ECG  and  relates  them  to  the  myocyte  membrane  currents. 


5.  An  In  Silico  Exploration  of  CN-Affected  Cardiac  Tissue 


The  ECG  modulation  caused  by  changes  in  the  myocyte’s  membrane  current  is  used  as  the  basis 
for  probing  CN’s  effect  on  cell  electrophysiology.  Effects  that  diminish  or  block  membrane 
currents  dictate  the  requirements  for  the  therapeutic  agent.  Computer  simulations  identify  the 
cellular  parameters  that  are  most  effective  in  controlling  these  changes. 

Available  models  (figure  2)  consider  only  indirectly,  if  at  all,  the  effect  of  metabolic  changes, 
especially  those  caused  by  disease,  on  cell  electrophysiology.  Myocardial  ischemia  causes 
anoxia,  acidosis,  and  an  increase  in  potassium  concentration.  In  some  respects,  it  is 
characterized  by  symptoms  analogous  to  CN  intoxication.  A  number  of  ischemia  models  have 
been  proposed  (Ferrero  et  ah,  1996;  Shaw  and  Rudy,  1997)  and  snapshots  have  been  given  of  the 
action  potentials  calculated  for  specific  levels  of  ATP,  based  on  experimental  observations 
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Figure  1.  Significant  CN-induced  changes  in  the  ECG  of  an  adult  male.* 


Reprinted  from  the  American  Heart  Journal  (Wexler  et  al.,  1947)  with  permission  from  Elsevier  Science:  New  York.  Additional  discussion  of  ECG  modulation  due  to  CN  is 
given  in  Katzman  and  Penney  (1993)  and  Salkowsky  and  Penney  (1995). 


Table  1.  Summary  of  significant  deviations  from  the  norm  found  in  the  ECG  of  a  CN-intoxicated  subject. 


Symptom 

Possible  Cause 

Primary  Membrane  Current  (I)  Affected 

P-wave  disappearance 

Auricular  arrest 

Depolarization  of  left  and  right  atria  by  INa 

T-wave  modulation 

Peaked  upright 

Inverted 

Polarity,  width 

Subendocardial  ischemia 
Subepicardial  ischemia 

Gap  junction  resistance 

Abnormal  repolarization,  IK 

ST-segment  elevation, 
shortening,  and  disappearance 

Ischemia,  hyperkalemia 

Ikatp  activation,  INa  inactivation,  Brugada 
syndrome 

QRS-complex  right-axis 
deviation 

Right  bundle  branch/left 
posterior  vascular  block 

Depolarization  currents  of  right  and  left 
ventricles 

Wenckebach  phenomenon 
(progressive  lengthening  of 
PR-interval  until  absence  of 
QRS-complex) 

Defective  conduction 
through  atrioventricular 
(AV)  node 

— 

Bradycardia  (FIR  <  50  BM) 

AV  block,  prolonged 
conduction  between  atria 
and  ventricles 

Decreased  ICa(L)  activity 

Tachycardia  (FIR  >100  BM) 

Increased  activity  of 
adrenergic  system 

If  in  sinus  node 

TdP,  VF 

Various,  including 
QT-interval  lengthening 

— 

Figure  2.  Effect  of  CN  on  the  cardiac  myocyte.  For  clarity,  not  all  dependencies  and  interactions  are  shown. 
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that,  in  anoxia,  the  ancillary  metabolic  changes  lead  to  a  decline  in  ATP  concentration.  Below 
an  ATP  concentration  threshold,  the  Ikatp  is  activated,  which  in  turn  affects  the  action  potential 
duration  and  shape  (Elliott  et  ah,  1989;  Light  et  ah,  2001;  Nichols  et  ah,  1991).  These  concepts 
can  be  extended  to  explain  the  ECG  changes  from  CN  intoxication. 

CN  inhibits  cytochrome  c  oxidase  activity,  which  in  turn  causes  a  decrease  in  the  mitochondrial 
respiration  and  exhibits  a  threshold  behavior  akin  to  some  of  the  mitochondrial  diseases.  The 
decline  in  ATP  concentration  opens  Katp  channels  (Borowitz  et  ah,  2001).  Also,  increases  in 
catecholamine  levels,  thought  to  contribute  to  the  genesis  of  arrhythmias  (Baskin  et  ah,  1987), 
intracellular  sodium  accumulation  (Kupriyanov  et  ah,  1996),  increases  in  inorganic  phosphate, 
lactate  production,  and  elevated  calcium  levels  (Kondo  et  al.,  1998),  are  noted.  This  report 
addresses  the  issue  of  ATP  concentration  change  effects  on  the  Ikatp  in  atrial  tissue. 

Atrial  Katp  currents  are  not  completely  understood  at  the  present  time.  Baron  et  al.  (2001) 
reported  that  pituitary  adenylate  cyclase-activating  polypeptide  (PACAP)  activates  the  Ikatp  in 
rat  atrial  myocytes  through  the  protein  kinase  A  and  protein  kinase  C  pathways.  PACAP  at 
1  nmol/L  activates  a  current  of  85  pA/pF  compared  to  26  pA/pF  with  5  mrnol/L  of  ATP.  A 
number  of  peptides  present  in  the  heart,  including  vasoactive  intestinal  peptide,  also  activate  the 
Ikatp,  though  not  as  effectively  as  PACAP.  Thus,  the  ultimate  model  will  have  to  take  several 
subsidiary,  and  heretofore  neglected,  events  into  account. 

Metabolic  impairment  that  activates  the  Ikatp  is  thought  to  have  cardioprotective  function  in 
ischemia  but  also  shortens  the  action  potential  duration  explained  on  the  basis  of  increased  [K+]0. 
Several  models  of  the  Ikatp  exist.  Ferrero  et  al.  (1996)  and  Shaw  and  Rudy  (1997)  have  included 
it  in  their  models  of  ventricle  electrophysiology.  These  models  propose  a  relationship  between 
the  current  density  and  ATP  concentration  and  give  snapshots  of  the  tissue  state  as  the  Ikatp  is 
activated.  Since  CN  intoxication  shares  a  number  of  characteristics  with  ischemia,  including 
anoxia  and  acidosis,  models  of  the  latter  are  convenient  starting  points  for  CN-affected  tissue 
simulation.  The  effect  of  the  CN  loading  on  the  atrial  tissue  is  modeled  by  expressing  the  Ikatp 
and  adjusting  the  potassium  and  sodium  concentrations  and  the  Ica(L)  • 

A  modified  atrial  model  (Nygren  et  al.,  1998)  is  used  to  demonstrate  that  modulation  of  the 
potassium  currents,  especially  when  the  Ikatp  is  present,  can  explain  the  AV  blockage  that  is 
expressed  by  the  P-wave  disappearing  from  the  ECG  during  CN  poisoning.  For  this  purpose,  an 
ATP-dependent  current  is  added  to  the  membrane  model  and  expressed  as  follows: 

Ikatp  =  facl  *  fac2  *  fac3  *  fATp  *  (Vm  -  EK),  (1) 

where  facl  is  the  channel  density  (~4.0  channcls/ptrT),  fac2  is  the  channel  conductance 
(-40.0  x  10  nS/cnr),  fac3  is  the  effect  of  the  presence  of  magnesium  ions,  fvip  is  a  function  of 
the  ATP  concentration,  Vm  is  the  voltage  in  the  membrane,  and  Ek  is  the  Nernst  potential. 
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In  the  resting  state,  the  ATP  concentration  is  taken  as  5  mrnol/L.  Three  cases  are  calculated: 

(1)  when  the  membrane  is  at  rest,  (2)  when  the  ATP  concentration  is  reduced  by  50%,  and 
(3)  when  all  of  the  ATP  has  been  used. 

For  the  calculations  reported,  the  monodomain  model  of  cardiac  tissue  is  assumed  to  be  valid. 
The  properties  of  the  cells  are  averaged  over  the  whole  domain  without  taking  the  anisotropy  of 
the  cardiac  tissue  into  account.  The  action  potential  propagation  is  detennined  by  solving  the 
following  equation: 


dv 

dt 


0,0 


1 


hon  0,0-  htim  0,0  +  -^  (^)V  o  (Dj  0)Vv0,  0) 

P  K  +  1 


(2) 


where  v  is  the  voltage,  i  is  the  current,  and  [>  and  k  characterize  the  tissue  with  D,-  as  the  diffusion 
coefficient. 

Obtaining  a  typical  solution  (i.e.,  the  time-resolved  electrophysiological  state  of  the  atria)  for 
300-ms  real  time  in  a  tissue  of  3  x  3  cm  required  7  hr  of  nondedicated  time  using  16  processors. 
Data  was  stored  at  5 -ms  intervals,  and  TIFF  files  were  constructed  using  voltages  as  functions  of 
position  in  the  two-dimensional  plane  of  the  tissue.  The  panels  displayed  were  constructed  using 
Tecplot,  version  9.01. 


6.  Results  and  Discussion 


CN’s  interference  with  atrial  tissue  energy  production  was  simulated  by  studying  the  effect  of  the 
changes  in  magnitude  of  the  Ikatp  on  the  action  potential,  thereby  understanding  the  processes 
leading  to  the  disappearance  of  the  P-wave  in  the  CN-affected  tissue.  Also,  ischemic  conditions 
were  used  as  the  hypothesized  state  of  the  tissue  with  the  Ica(L)  partially  blocked. 

Figure  3  shows  the  atrial  action  potential  at  several  concentrations  of  ATP,  starting  with  the 
nonnal,  homeostatic  and  corresponding  to  [ATP]  at  >6.0,  at  2.0,  and  at  1.0  rnmol/L,  respectively. 
Notable  shortening  of  the  action  potential  is  observed  from  300  ms  down  to  10  ms  at  [ATP]  = 

1 .0  mmol/L.  In  tandem,  the  slope  of  the  repolarization  changes.  The  maximum  value  of  the 
membrane  voltage  drops  from  ~20  mV  to  the  single-digit  values  under  conditions  of  increased 
[Na+]i  and  [K+]0,  that  is,  accumulated  sodium  ions  inside  of  the  cell  and  potassium  ions  outside  of 
the  cell.  When  the  potassium  and  sodium  are  held  at  homeostatic  values,  the  action  potential 
does  not  reach  the  threshold  voltage,  indicating  impending  shutdown  of  the  cell. 

Drastic  shortening  of  the  action  potential  cycle  is  observed  under  conditions  of  elevated 
potassium  and  sodium  concentrations  (figure  4),  with  a  slightly  elevated  magnesium 
concentration  at  3.5  mmol/L  (approximating  CN-toxicity  conditions)  and  a  75%  block  of  the 

ICa(L)- 
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Figure  3.  Effect  of  ATP  concentration  changes  on  the  atrial  action  potential.  The  top  curve 
denotes  the  homeostatic  condition.  The  middle  curve  gives  the  condition  for  [ATP] 
=  2.0  mmol/L,  and  the  bottom  curve  gives  the  condition  for  [ATP]  =1.0  mmol/L. 


Again,  the  homeostatic  cycle  length  (CL)  of  -300  ms  is  shortened  to  10  ms  at  [ATP] 

=  2.0  mmol/L,  with  CL  =  <10  ms  at  [ATP]  =1.0  mmol/L.  The  action  potential  maximum 
voltage  value  stays  positive  though  considerably  lower,  -8.0  mV,  than  the  normal  value.  Also, 
the  voltage  never  returns  to  the  baseline  value  of -74.0  mV  but  approaches  -68.0  mV. 

At  elevated  [Na+]j  values  of  17.0  mmol/L,  the  action  potential  hyperpolarizes.  The 
hyperpolarization  disappears  as  the  [K+]0  is  elevated.  Under  these  conditions,  the  normal  resting 
voltage  is  never  reached. 

Figure  5  contrasts  the  normal  action  potential  with  the  case  when  [ATP]  =  2.0  mmol/L  and 
[Na+]j  =  4.5  mmol/L,  a  subnormal  condition  when  [K+]0  is  held  at  4.5  mmol/L.  Again,  the  CL  is 
shortened,  the  maximum  amplitude  is  lowered,  and  the  slope  of  the  repolarization  changes 
drastically. 

The  reduction  in  time  available  for  the  depolarization  to  take  place  (figure  6)  as  well  as  the 
voltage  amplitude  observed  under  these  conditions  makes  a  strong  case  for  the  diminution  and 
eventual  complete  disappearance  of  the  P-wave  from  the  ECG.  Also,  blockage  of  the 
depolarization  wave  is  a  plausible  scenario. 

The  study  described  in  this  report  addresses  only  indirectly  the  effects  of  calcium  ion  flow  and 
concentration  changes  by  including  a  partial  block  of  the  Ica(L)-  The  rundown  of  the  ATP  energy 
supply  for  the  sarcoplasmic  reticulum  pump  that  transfers  calcium  ions  into  the  extracellular 
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Figure  4.  Atrial  action  potentials  under  ischemia-like  conditions,  caused  by  CN  toxicity,  that 
increase  the  [Na+];and  [K+]0  levels  in  the  tissue.  The  homeostatic  condition  (top)  is 
compared  with  conditions  of  [ATP]  =1.5  mmol/L  (middle)  and  1.0  mmol/L  (bottom) 
at  [Na+]|  and  [K+]0  of  17.0  and  12.0  mmol/L,  respectively,  with  a  slightly  elevated 
[Mg2+]  of  3.5  mmol/L. 


Figure  5.  Contrast  of  a  normal  atrial  action  potential  (top)  with  one  when  [ATP]  is  depleted  to 
2.0  mmol/L  and  both  [K+]0  and  [Na+];  are  at  a  subnormal  level  of  4.5  mmol/L  each 
(bottom). 
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Figure  6.  Atrial  tissue  with  [K+]0  elevated  to  12.0  mmol/L  ([top]  initial  resting  voltage  is  no 
longer  reached)  and  [Na+];  elevated  to  17.0  mmol/L  ([bottom]  the  resting  voltage  is 
depressed,  but  the  trace  remains  above  the  original  value). 


space  against  a  concentration  gradient  leads  to  calcium  overload.  This  may  lead  to  abnormal 
calcium-induced  currents,  which  are  precursors  to  arrhythmia.  Changes  in  action  potential 
morphology  express,  and  possibly  contribute  to,  the  modulation  of  the  Ica(L>  The  threshold  for 
the  activation  of  this  current  in  myocytes  is  around  -25  mV.  Half-maximum  deactivation  is 
observed  at  or  below  -20  mV.  One  of  CN’s  effects  on  action  potentials  is  the  narrowing  of  the 
area  under  the  curve  with  the  decline  in  ATP  availability.  This  is  signaled  by  an  increase  in  the 
angle,  measured  from  the  horizontal,  from  the  maximum  of  the  voltage-vs.-time  plots.  The 
narrowing  indicates  a  decrease  in  the  time  that  the  calcium  ion  channels  remain  open  and 
operative  and,  thus,  the  amount  of  calcium  ions  that  are  transferred.  Consequently,  the  total  ion 
flow  is  reduced. 

The  electrophysiological  changes  caused  by  CN  intoxication  of  the  atria  lead  to  profound 
changes  in  the  membrane  currents  needed  to  maintain  the  P-wave  in  an  ECG.  Indeed,  the  effect 
of  the  CN  is  to  render  the  membrane  current  levels  insufficient  for  maintaining  the  required  wave 
shape  and  duration.  These  results  suggest  that  pharmacological  intervention  could  seek  the 
restoration  of  ion  concentration  levels  as  a  possible  approach  to  restoring  homeostasis. 

The  simulation  of  cardiac  toxicity  with  ancillary  markers  generated  by  the  pathology  requires 
consideration  of  a  tissue  domain  large  enough  to  allow  the  cellular  processes  to  establish  a  basis 
that  leads  to  wave  propagation  and  other  manifestations  of  the  modulation  of  the 
electrophysiologic  processes.  A  myocyte  is  typically  100  x  20  pm,  thus  even  a  small  patch  of 
tissue  consists  of  a  very  large  number  of  cells.  A  simulation  based  on  cellular  processes  requires 
sizeable  computational  assets. 
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The  computational  resources  required  for  the  calculation  of  a  1  x  1  cm  piece  of  tissue  is  shown  in 
figure.  7.  When  fewer  than  16  processors  are  used,  the  computational  times  become  prohibitive, 
jumping  from  1830  s  for  16  processors  to  3356  s  for  8  processors  and  6566  s  for  4  processors. 
The  scalability  of  the  algorithm  is  apparent.  The  processing  time  with  64  processors  was  611s, 
underlining  the  advantage  of  parallel  processing. 


Figure  7.  Run  time  as  a  function  of  number  of  processing  nodes  on  the  SP-4. 
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